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I . INTRODUCTION 



A. BACKGROUND 

This paper describes the concept and construction 
of a hard wired analog device, pictured below, which 
simulates the operation of a gas turbine engine. 




Recent work on gas turbine simulation has focused 
largely on digital or hybrid computing. To accommodate 
the nonlinear elements of a gas turbine simulation in this 
fashion requires a perturbation or piecewise linear 
technique. A pure analog approach offers a continuous 
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solution to these nonlinear characteristics, which facili- 
tates a simulation in real time. 

At this time, analog computation of complex non- 
linear systems is not feasible on devices such as the EAI 
680 computer currently installed in the M.I.T. Joint 
Computer Facility. With the miniaturization of both digital 
and linear functions into the now-common integrated circuit 
or IC chip, it becomes possible to design extremely compact 
analog devices. The IC component cost of such a simulation 
is rapidly approaching the cost of purely passive electronic 
elements. 

The modeling approach to the analog simulation, 
physical modeling, follows the work of Markunas^ closely. 
Physical modeling attempts to explain the engine char- 
acteristics within a framework of physical princples and 
laws which describe the functions occurring in the real 
machine. Simplicity is generally a strong point of the 
model . 

The model which has been implemented in the analog 
simulation represents a general machine "invented" by 
Markunas, who estimated the necessary characteristic 
relations between the independent variables. Markunas 
conducted digital and hybrid simulations of his "general" 
gas turbine, but a full analog implementation was not 
investigated. 
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This paper is essentially a design manual which 
traces the development of a hard-wired analog simulator. 
It is intended to give the reader a brief overview of the 
analog devices currently available, and their application 
to a computational circuit. To facilitate future work 
with the simulator, a detailed user's manual has been 
written. 

This document contains explicit details relating 
to the construction and operation of the simulator. The 
user's manual, the simulator, and associated hardware 
are under the supervisor of Professor Henry M. Paynter, 
Department of Mechanical Engineering. 

B. DESIGN CONCEPT 

The following characteristics were desired in the 
hard-wired simulator: 

1. Real time, non-linear third order simulation 
of the dynamics of a single shaft gas 
turbine engine, from start-up to overload 
conditions . 

2. Cost to be under $600. 

3. Size of the device such as to be readily 
portable and convenient to set-up and operate. 
Layout and operation of the simulator 
to be comparable with typical industrial 



4 . 
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and military engine control consoles. 

5. Accuracy of the steady-state simulation to 
be considered of secondary importance 
compared with the representation of large 
scale transient effects. 

6. The ability to investigate various fuel 
control schemes in a closed loop control, 
and the ability to operate with different 
torque-speed load characteristics. 

7. Utilization of the general purpose, low 
cost IC operational amplifier in a hard 
wired configuration to accomplish all 
mathematical operations. 



II. MODELING 



A. GENERAL APPROACH 

The gas turbine model selected by Markunas will be 
briefly explained. The modeling approach is broken down 
into: steady state and dimensional considerations, 

frequency analysis, and dynamic equations. 

Markunas' estimated performance parameters and 
the dynamic equations are combined in state-variable 
format. These equations appear in Appendix II and form 
the basis for the analog signal flow used in the design. 

B. STEADY STATE AND DIMENSIONAL ANALYSIS 

3 

Horlock establishes non-dimensional functions in 
describing the performance of compressors and turbines. 

A set of relations are developed to describe the 
performance of a turbomachine in its simplest terms. 

In the subsequent development, gas turbine state 
points are used, numbered in the following convention, with 
total values indicated by a zero in front of the state 
point: 
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The following symbols are used for dimensional 
analysis in the compressor: 



Symbol 


Compressor 


Dimensions 


fhrz 


Inlet stagnation pressure 


n/LT^ 




Inlet stagnation density 


r i/l* 




Outlet stagnation pressure 


/ lt' z 


To*Z 


Inlet stagnation temperature 






Efficiency 


— 


u> 


A characteristic linear dimension 




Mo 


Stagnation temperature rise through 
compressor 




M 


Rotational speed 


«/T 


VVk 


Mass flow rate of gas 


bA /t 
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Symbol 


Compressor 




Dimensions 


jn 


Molecular wt of gas 




— 




Absolute viscosity of gas 




m/lt 


cp 


Specific heat at constant 


pressure 


l 7 /r‘‘e> 


'zr 


Ratio of specific heats 




— 


R 


Gas constant 






Q.O.-Z 


Speed of sound of entering 


gas 


L./'T' 


'll 


Torque input 




ri C'/t^ 


In the 


compressor, if lx 


and /\To are 


> chosen as 



variables of interest, one can write 

'Fcyss ) ) Tl_ to ■=.-£’( Fci-z. 1 Tg>-z. ) > fP, rrv ) ^ a -2 > ^ ) > Q c> 

Assuming a perfect gas, fixing the gas constant , G2., choosing 

, and "Tc>x thus constrains ^ 03 . (’oz can be expressed in 
terms of (5 ) ; the original relation becomes 

CvTi ^ ,N,G>,<r?v 

Taking non-dimensional groups, 

'r: £>To 






^ G 



h*P m J rTc » 2 ^ P02D 
S/tSi > R. z D" 2 ^'SxSz 



vhz , rc-z 

With constant values, and ^ can be eliminated: 

✓ 

P<S~3 ^7 XL_ __ . &To ~ c ( NP . . I'ozO 

iS _ 1 v ' iJ? ...r-3 ’ “ V > — X \ r 5 — — 



V-fe'at 



Toz 



'friz > D' 






Multiplying the last two non-dimensional parameters to 
create a new parameter. 
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FWi , £jo --P/VlD i rr- . 

' R*D* ’tT^ to. D* D 



It can be seen that the third parameter, is a 

"flow Reynolds number". 

Assuming that the compressor operates above some 
critical Reynolds number, and that the dimensions of a 
particular machine are fixed, the equations reduce to: 



Pb>^ ^ ) TZ ^ - -f f H ) 

Pci'Z. ^ ir^-2L TVz, \ V 4 Tg' 2L. ^ Pcvz_ / 

In order to present meaningful performance relations, 
compressor data should be corrected following the above 
parameters . 

4 . 

Identically for a turbine, it can be shown that 
the same relations hold and that the parameters of 
interest in turbine performance"^ are d*V> ) sv\ ) I 'Z. ) • 

All inlet conditions assume values at state point 4, and 
all exit values are at state point 5. Torque across the 
turbine is ‘"Z . As shown before. 




ATo _ £ 





The characteristics for a given turbomachine can 
now be specified as functions of dimensionless or pseudo- 
dimen s ion le s s group s : 
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Compressor : 

^PjL 

Fo-2 

Turbine : 

Fc^\ 



Xo^ 

To* 




■£ 



R>^ \ 

^ST-i 4 TcTSL / 



„ / ^5 ..Nt 

^ k pe> -l , 'frS^ 




) 



< 



Po-r \ 

poq > \T'fcA / 

ft*3T JfcsL \ 

Pb4 ^ \fTvi ) 
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C. FREQUENCY ANALYSIS 

Modeling of the processes which take place in a gas 
turbine can generally be classified as mechanical, fluid, 
or thermal in nature. To reduce the scope of the problem 
to a more manageable field, only the most significant 
energy domains and dynamics were considered. Markunas 
defines three dynamic regimes for the gas turbine: 

1. Low frequency, (f ^ .1 hz) where rotary 
inertias and thermal capacitances dominate 
the dynamics. 

2. Medium frequency (f ££ l hz) where rotary 
inertias and fluid capacitances dominate. 

3. High frequency (f^t 100 hz) where the fluid 
dynamics are dominant. 

Markunas established that rotor dynamics and fluid 
capacitances were the primary influences on gas turbine 
dynamics. Thermal capacitances and fluid inertias were 
neglected. 

D. DYNAMIC EQUATIONS 
1. Assumptions 

The gas turbine model is assumed to operate with a 
perfect gas in one-dimensional flow. Components are 
treated in lumped parameter fashion. 



-17- 



2. Rotor Inertia 



For 


the 


rotor, a summation of torques give the 


following: 


Nier 


— 


3=0 
TT l45 


where 


Mcx 


= 


angular acceleration, rpm/sec 




= 


2 

rotary inertia, ft-lbf-sec 




= 


load, compressor, or turbine torque, ft-lbf 


3. 


Fluid Capacitance 


For 


the 


one-dimensional lumped parameter fluid, the 


following 


symbology is used : 


e 


= 


fluid density, lbm/ft 


* 


= 


axial distance along duct, ft 


v* 


= 


mean axial velocity of fluid, ft/sec 




= 


ratio of specific heats = Cp/Cs/ 


Cp 


= 


specific heat at constant pressure 


Cv 


= 


specific heat at constant volume 


N\ 


= 


Mach number = VV / 


Cv 


= 


Local sonic velocity = \l ^ 


<j 


= 


internal energy of the fluid in ft-lbf/lbm 


K 


= 


enthalpy of the fluid in ft-lbf/lbm 


y<*L- 


= 


lumped fluid volume 
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A general mass conservation equation for the one 
dimensional fluid stream may be written: 

Integrating over a constant area duct, with fluid pro- 
perties constant at a given axial distance or dimension: 

2b Jr ^ (VA'v'*'') 



“ 2 >^ 



Integrating over the length of the duct 



by definition. 









(f = o 






— o 



is approximated by an average density, 
^ • Lumped fluid volume 



since 



4 (t) - 
* 3 + 



\/ot 



c 



\TAx - 
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where the subscipts 1 and 2 identify inlet and outlet 
conditions across the engine, respectively. 

Writing a similar equation for the conservation of 



momentum neglecting viscous shear and body forces. 







-V 



/f Vy? 




= Cs 



Integrating over area, 



, 



(&&■') * 










o 













- c> 



Integrating over length of the duct 
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i- 

7 



Since total pressure. 



IS— 



- ?(l y ) v ' 1 * 



d_(rn") = 
dr 




Fo 

JT7*± 



-L 



2 . 



Assuming V = 1.4, and M £ .5, the momentum equation can 
be written: 



^ ("Po.-Po^ 



Similarly for energy conservation with no external 



work or heat addition, 

- z 






_ 













— o 



Integrating over area and length, 

U 






fACuvVy *) 

-I 



r *!^>L 

■ze> c ' 



C5 
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JL 



( V<0 ^ 



( 2(0 vvV \ 



— m 



where 



Since 



Ko H 



u + vk 

‘2*0 



*2 

C 



"^ < 5 C - 

Substituting , 






O 4- V / x.' Z -- CvT" r _L (V~ 0 ^ *N 

-2C^ C v 7 



M? CvU ) 



«■ 

\ JL+- 



l -4- Jn/^‘ 



- — /wJoi.-«^VzTc 



vcu\ 



oz 



where 



la = T 4- (X f ^ M 2 '") 

With ~ *3" the result can be simplified to: 

*(f r * S ) ' tL 




E. STATE EQUATIONS 

The three conservation equations are now: 
1. Mass 



^(f') - 4_ ^ 
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2 . Momentum 

3 . Energy 

Hr 

Using the definition of a polytropic process, 

_ n 

■= const ( it follows that 

r> 

(coK^r) 

T 0 - Ccota^fe " = (Tc^sr^fo 

In order to get the mass and energy equations into more 
convenient form, 




( C£>Nt < £T' S ) ( h ^ 







I 
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With the result that the equations now be more: 




Where 1 and 2 represent inlet and outlet across the 
duct. Because of lumped parameter considerations, machine 
model parameters can be substituted as follows: 



One Dimensional 

O 

AV 

crs-s 

TSa. 

Tc^ 

Po 



Machine Model 

<V\^, 

Tc^ 

To-s 



Substituting notation for the first and third 
differential equations. 
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Tc>4- 



^To4- 
^ Fo- 5 Vot. 




Introducing the torque equation as the third system 
equation , 




where 

= fuel lower heating value, Btu/lbm 

% 

= fuel flow, lbm/sec 
» ^ = load torque, ft-lbf 

kMO are chosen as variables of 

convenience. 

For an isentropic process, n = 1.4. Markunas felt 
that n-l/n should be replaced by 1.0 in the temperature 
equation. 

F. ENGINE DATA, ASSUMPTIONS, AND CHARACTERISTICS 

1. Introduction 

In Markunas' analysis a single shaft machine was 
invented from a number of sources. Several real world 
machines exist which have design parameters close to 
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Markunas' s engine. The recent Garrett 
its locked shaft mode is a reasonable 
2. Design Point Specifications 


GTPF 990 engine 
example . 


Pressure Ratio 


5.0 


Compressor Isentropic Efficiency 


1280°F 


Turbine Isentropic Efficiency 


88% 


Compressor Discharge Mass Flow 


100 lbm/sec 


Inlet Pressure 


14.175 psia 


Inlet Temperature 


80°F 


Rotor Speed 


7200 rpm 


Output 


6948 HP 



Mar kuna s plotted characteristic curves for his 
compressor and turbine, which appear in Appendix I. 
Markunas then mathematically "fit" each characteristic 
curve. These approximations follow the functional 
relations established for non-dimensional and pseudo 
dimensionless groups. 

3. Compressor and Turbine Characteristics 
Markunas defined the following variables and char- 
acteristics for his compressor: 



PR.C 



pressure ratio 
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TRC - 

xwc - 

XNC 

XMC • 

NC - 

ec ^ 
cc = 
oc - 

EC *■ 
PC - 



To^ /To-z 



temperature ratio 




scaled compressor 
mass flow 



scaled compressor 
speed 



vfrS 



07 



( yv >%\|ri 



I?' 2 ') 

P<sTZ CX'Sl < 54*4 



.~7S C*M<0 * .^(vMc) 



4 . 

"^/O V -^.Ci 0**^) 



-s-% 

. *>l£> C~z — XMQ*) 

s[7 - CXMC — /O^ 

4-1 

±.,o +■ XNC t “S.fcCxNK') 

4,0 ur ppC £ gc 

- 2 . ,0 tF P(2C *> £ C 
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CiEC - 


*y — -V- 
' ^ 


C 


y>wc = 


A£ 


[l-(pcc-i.n] 

L - Be J 


EFFC = 


Dcj 


1 -^ec -fgscN 

- N EC -JL > 


TRC ^ 


1 — ' ^rEC 

X+KPRC') -L 



EFEC 



Similarly, for the turbine 



Per - 



rRr 

ywr 

V.MT 

XMT 

ATT 



X^/Toa 



s/r c 



VYVV-V Va 4 



/ 






V ^ 



vJTcy 












Po>^g / s /^4 ) 
vi Tcyv ^ oe%v<-^ 



Mc.\ 



^4 / \^4 



1 . 



~S-L 



( 



X 

sKrt-^y 
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cr - 


n 

r 


4 — | 




or 




JL- (jL- *MT ) 

mm m 




ET 


1 . 






<=-ev -= 






V.VNT 


Cs.t[ 


a - C PR-rf T 1 CT 


EFFT = 


I 

X - 




-2 


TW 


X - 


(_ M 1_ - PRC _ 

• EFFT jj. _ CmT 


c^t 

•) 



Additional Machine Assumptions 

1. =r (no duct loss) 

~ ^ (no combustor losses) 

P03 - FM 

The machine characteristics are combined with the 



differential equations to produce a set of system state 
equations. The expanded equations appear in Appendix II. 
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III. THE OPERATIONAL AMPLIFIER 

A. INTRODUCTION 

The fundamentals of the operational amplifier, (OA) 
are briefly presented to establish the necessary analog 
"tools" to proceed to a simulator design. 

B. BASIC OPERATION 

An ideal OP amp can be modeled by the following 




with no differential voltage between positive W ) and 
negative ( V ^inputs, both emitters are equally biased, and 
collector currents Q±. and C^ 2 are equal. The differential 
voltage output between 1 and 2 is also zero. Driving V 
more positive with a signal will cause the current through 
<^.to increase. The current through must decrease by a 
proportional amount, so that in the differential voltage 
output between 1 and 2, 1 gets driven more negative with 
respect to its initial voltage, and 2 becomes more positive. 
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torns for a small signal applied to one of the input terminals, 
a high voltage gain and current "buffering" takes place as 
a result of the transistor characteristics. 

With the ideal OA, the input current is essentially 
zero, and the voltage gain is assumed infinite. 

C. THE INVERTING OA 

Looking at the following diagram, 




It is assumed that the impedance of the amplifier 
, therefore Tc% c> - 

Closed Loop Gain, from a summation of currents: 






C Vo — Vo. ^ / ' 



xx - Cvr - s/O /'Ei • i-o. ~xto 

- Vo / S' 3 t Vn/ V i /Wx — ( zrx 



V 

\/r 



- Qw '2-a/C^j. 
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If a, the amplifier gain, is large, 

N/c, - jj ' — -z 

V,' ~ 

It can therefore be seen that the inverting OA inverts 
and provides a gain proportional to the input and the 
signal feedback resistances. 

D. INVERTING SUMMER 

/ 

By joining multiple input currents to the input side 
of the OA, an inverting summer can be constructed with 
variable gains on each input. 




By current summation. 
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E. INTEGRATOR 



Adding a capacitor to the feed back loop 
adds the following constitutive relation to the 
loop : 




of the OA 
feedback 



Looking at a typical circuit. 




X* -X\=> 



£sV? _ Xk2 
t± f c 



zWo ^ 


- v; 








•+* 

y 


Vo ~ 


— 1 . 


( v/; dr 




RiC 


\ 
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F. LOGARITHMIC GENERATOR 

Using a diode in the feedback loop enables the 
logarithmic nature of the diode voltage - current relation 
to be included in the loop. The diode characteristic 
can be modeled by the following equation: 




where : 



T.o 

Vo 

X 

V 



junction saturation current 
junction saturation voltage 
bias current 
bias voltage 



A typical characteristic appears as: 



tv\|\ 



To 

* 



v 
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It is important to note that the saturation voltage is 

directly proportional to the absolute temperature. 

A typical log circuit can be constructed as follows 




A diode - connected transistor is used to perform 
the diode function, as transistor emitter-base junctions 
have more dynamic current range than diodes in general. 

By the following diagram, an input V»' across Rl causes 
a current flow through the " transdiode" . The bias voltage 
appearing as varies as the log of V*’ . 



CoeC.feMt 



Tc,. 

Tfe 













-V— 

Vo 



v; 
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G. ANTILOG GENERATOR 

In a similar fashion, input voltages or currents may 
be exponentiated in the following circuit : 
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IV. ANALOG COMPUTATION 



A. INTRODUCTION 

To implement the non-linear system equations the 
following functional operations are performed: addition, 

subtraction, integration , multiplication, division, 
exponentiation to a fixed power, and exponentiation to 
a variable power. The first three analog operations have 
been briefly discussed. To implement the higher functions 
the use of logarithmic circuits must be made. 

B. THE PRACTICAL LOG CIRCUIT 

There are at least three major considerations in 
the construction of a practical log or antilog circuit: 
cost, complexity, and accuracy. Three different 
architectures are considered to accomplish the log function 

The first choice is the previously discussed simple 
log/antilog circuit. Using low quality components, cost 
of the log circuit is under one dollar. Thermal stability 
is considered poor. A test circuit and transdiode 
characteristics are shown in Appendix III. Test data 
indicated that the accuracy and thermal properties of 
this circuit were marginal. 

A second choice is a compensated log circuit as 
diagrammed below, which offers increased thermal stability 
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In the diagram, the collector of Q 2 . is kept at virtual 
ground by the OP amp, and the collector current through 
is equal to a constant reference current provided 
by • It can be shown that the output voltage is of 

the form Vo “=? ~ C 3r ~L ^ where the thermistor 

j — uL / 

provides temperature compensation. One of the drawbacks 
to this scheme is that the number of components necessary 
to generate the log function has gone from 3 to at least 
9, with proportional space and cost increases. This 
configuration was not evaluated. 

The third and most expensive choice is to utilize 
a dedicated logarithmic device. The particular unit 
tested was a Texas Instrument log amplifier. This is a 
dual log generator on a single DIP package, at a cost of 
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approximately $2.50. With appropriate circuit manipulation, 
the desired log/anti-log operations can be performed. The 
TI log amplifier was evaluated in depth. Results are 
documented in Appendix III. The log amplifier was the 
foundation of the simulator design. Other dedicated log/ 
antilog modules exist which offer increased dynamic range, 
temperature stability, and accuracy, but alternative 
devices were an order of magnitude more costly than the 
TI log amplifier. 



C. 



EXPONENTIATION - ALL CASES 



In order to achieve either fixed or variable power 
exponentiation, the use of logarithms is basic. The 
signal flow for a fixed power exponent (power function 
generator or PFG) circuit is diagrammed below: 
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By a judicious selection of values for and V^2’ 
current can be multiplied or divided using an operational 
amplifier in the ratio: 



To be able to exponentiate to a variable power 
requires a cascaded log/antilog arrangement. The use of 
a segment of the state equations is made to demonstrate 
the concept. It is desired to generate the following: 

c<: 

wwc ■= vc ( l ~ C p rc -c ) ^ 

B.C 

where 1C represents the scaling of several fixed parameters, 
and \W C ) > 0C and ^ are all problem variables. 

The circuit selected must process: 



L_C*-» LOGt 
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As a result of thermal instabilities, the interior section 
of the signal flow was replaced with a four-quadrant 
multiplier. 

The multiplier essentially compresses the separate 



functions into a single chip. By circuit manipulation, 
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a multiplier can also divide, square, and take square 
roots. A range of multiplier devices currently exist, 
with the price of the least expensive about $3. 

Test data on the multiplier used in the simulator 
was collected in Appendix III. To allow as compact a 
circuit as possible, the multiplier was frequently 
utilized in its dividing and exponentiating modes. 
Complete device specifications are included in 
Appendix VI . 
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V. PRELIMINARY DESIGN 



A. INTRODUCTION 

The simulator design proceeded in a parallel fashion. 
In an iterative process, estimations were generated for the 
electronic layout, which allowed the cabinet and control 
panel to be sized and designed. 

The simulator is configured to resemble a real-world 
engine control panel. Additional concepts for the per- 
ipheral design which were not incorporated were: 

1. Proportional fuel control schemes with 
variable gains. 

2. Torque-speed loads characteristic of 
pumps, propellers, generators, and other 
real world devices. 

3. Overspeed alarms and stall audio noise 
circuits . 

B. ELECTRONICS 

1. Logic Circuit 

The logic processing of the electronic circuitry 
was broken down into 3 separate grid structures, which 
appear as design drawings in Appendix V. Use was made 
of bread-board panels to avoid the inflexibility of other 
techniques. Concurrent with testing various mathematical 
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operations, the circuit design was modified to balance the 
competing parameters of simplicity, size, accuracy, and 
cost. The appropriate equation variables were normalized 
about their design point values, and then scaled to 
provide a circuit signal that was approximately 75% of 
the maximum rates input for the device. The steady 
state design point values for all variables were established 
by a DYSYS program and used for the normalization. 

Additional notes on the wiring techniques and breadboard 
layout are included in Appendix VI. 

2 . DYSYS Simulation and Scaling 

A DYSYS simulation of the system equations was run 
as it appears in Appendix IV. Significant comments on 
the program are : 

1. The conditional value of FC was changed 
to a fixed value of 3.0 to facilitate the 
circuit logic. 

2. Markunas estimated a parabolic "pump-type" 
load torque, which appears as . 0000977 (3) **2 . 

3. The modeling parameter DT is called TT in 
the program to avoid confusion with the 
DYSYS time-step variables, DT. 
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The scaling procedure started with the following 
listings : 



a. Constants 
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c. 



Input Variables at Design Point 





±, L&CT7 

-5*o i -"» - u^v- ' 



d. State Variables at Design Point 
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The DYSYS simulation was run to verify Markunas's 
simulations. Use was made of the program to assist in 
adding voltage and current values to the circuit design, 
and to evaluate the effect of fixing the discontinuous 
constant, FC. FC was conditionally assigned the value 
of 2 or 4 in Markunas work. The steady state effect of 
assigning FC as 3 was investigated and adopted. 

The DYSYS simulation was also run off - design 
to generate a test case for the logic circuit. The 
simulation was run at approximately 50% fuel flow and 
steady-state values were generated for all elements of 
the equations. 
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Some difficulty was encountered in determining the 
correct value for the lumped parameter fluid volume. The 
value which appeared to be consistent with Markunas' data 
was 144 ft^. An approximation based on real world data 
gave 15 ft^ as a more likely figure. Since the volume 
represents a time constant scaling for the fluid capacitance, 
it has no significance for the DYSYS steady state values. 
Changes in the volume assumed can be readily investigated 
in the simulator dynamic response. 

The scaled voltage values and the system variables 
they represent are drawn on the 3 circuit plans in 
Appendix V. 

C. PERIPHERALS 

The main elements of the remaining design work 
included the chasis, power supply and distribution, 
instrumentation and control hardware, and breadboard 
architecture. The design drawings and sketches included 
in Appendix V illustrate the final peripheral design. 

The photograph below shows the internal layout of the 



simulator. 
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Appendix VI, specifications, contains additional 
data on the peripheral equipment, its wiring, and operating 
parameters. 
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VI. IMPLEMENTATION AND CHECKOUT 

A. PROCEDURE 

To establish the characteristics of the various 
devices and to investigate the methodology of the design, 

10 tests were conducted. The significant results are 
discussed here in the context of the design process. All 
test data is included in Appendix III. 

B. LOG GENERATION 

Tests 1 and 2 determined the voltage - current 
relations for the transdiode - connected OP amp, using 
both 2N2222 and 2N2904 transistors. It became apparent 
that the allowable voltage values of the log circuit 
were so narrow that only current variations could be used 
for the signal, over about a 1 1/2 decade range of values. 

C. POWER EXPONENTIATION 

Tests 3, 4, 5, and 6 were designed to varify the 
fixed power exponentiation circuitry. A workable 
configuration was wired in test 6, using a current - 
dividing resistive network. 

Tests 7, 8, 9, and 10 were conducted to validate 
the variable - power exponentiation procedure. The inability 
of the transdiode - connected OP amp to provide a stable 
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cascaded log circuit led to the investigation of a 
4 quadrant multiplier in Test 8. 

The Exar 2208 multiplier characteristics were 
investigated to generate hookup diagrams, trimming 
values, and allowable voltage values. 

With Test 9, a simple log/multiplier/simple 
antilog circuit was wired and observed. Stability was 
significantly improved, but was still considered below 
the accuracy level desired. 

In Test 10, a Texas Instrument 441 log amplifier 
was tested to determine its log/antilog characteristics. 

The log- amp was incorporated into the variable expon- 
entiation circuit and the stability evaluated. The 
steady state accuracy of the log amplifier in the 
variable exponentiation circuit was somewhat better than 
the simple log configuration. Since the dedicated 
log amplifier offered better compensation characteristics 
than the simple log configuration, it was substituted 
at all points in the circuit design. 

D. GENERAL CIRCUIT DESIGN 

Concurrent with the test sequence, the overall 
circuit design was changing as new test data was generated. 
The more significant problems which caused iteration of 
the circuit are briefly discussed. One immediate limitation 
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exists using the TI log amplifier module. Looking at 
"TL 441 Transfer Characteristics" included in Appendix III, 
the log amplifier only has a useable output range of 0 - 
450 millivolts. This can be doubled by wiring an OP amp 
to the input side of the log amplifier to use the dual 
input capability of the device. Still, the output accuracy 
is limited, since various offsets for multipliers and OP 
amps can easily reach 10 millivolts. 

As a general rule, it was decided to trin the 
circuit at a limited number of points, rather than dealing 
with the offsets required for each component. 

It was intended in the final circuit design to 
operate all components at one dual-polarity supply voltage. 

This choice nominally would be + 8 volts, the maximum 
rating of the most sensitive device, the TI log amp. Test 
10 demonstrated that significant nonlinearities were intro- 
duced into the multiplier characteristics if the supply voltage 
were reduced. Three solutions were postulated: 

1. Run the multipliers at + 15 volts, and the 
OP and log amps at + 8 volts. 

2. Compensate the multiplier to have linear 
characteristics at + 8 volts supply. 

Replace the multiplier modules with log 
amplifiers to process the same function. 



3 . 
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The first choice was selected in an effort to keep 
the circuit design as simple as possible, in as much as 
the single - voltage supply was more of a packaging 
consideration than an electronic one. 



V 
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VII. CONCLUSIONS AND RECOMMENDATIONS 

A. FEASIBILITY OF THE DESIGN 

From the test data collected, a real time simulation 
can be conducted with the logic circuits wired as designed. 
The techniques developed for the implementation of non- 
linear differential equations should be applicable to a 
wide range of system models. 

Markunas defined his engine in essentially non 
dimensional or pseudo-dimensionless groups. It was not 
possible to investigate the ability of the simulator to 
represent different geometric machines, however, it is 
felt that future work will verify the utility of the 
simulator and the nonlinear model. 

B. DESIGN REFINEMENTS 

The simulator as designed could be significantly 
enhanced with more general torque and fuel control 
schemes. If the simulator could demonstate sufficient 
generality in modeling the behavior of real world gas 
turbines, an autonomous "observer" circuit could be 
configured for a wide range of engines. This observer 
device, using the logic developed in the simulator, could 
be an extremely compact unit. It could provide improved 
control of fuel, speed, temperature, and pressure 
dynamics. Such a circuit could be incorporated in one 
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or two monolithic LSI chips, to be externally trimmed for 
a particular engine or operating condition. 
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APPENDIX I 

MACHINE CHARACTERISTICS 



AS 



- 56 - 







- 57 - 







±/S" 



- 58 - 



TT 




SCALE P C*MC) 



SCM-£C> 

_ 



- 59 - 




OeO 



- 60 - 




0 

"I 






- 61 - 




'SC/M-ED -spcgo Cxnt^) 



- 62 - 



APPENDIX II 
STATE EQUATIONS 
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APPENDIX III 
Test Data 



A. TEST 1 AND 2 

19 

From reference data, the transdiode voltage current 
characteristics was estimated for various values of 
saturation current. The log circuits diagrammed subsequently 
for the 2N2222 and 2N2904 transistors demonstrated low 
base-emitter saturation currents, thus making the exponential 
characteristic curve more abrupt. 
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Test Data 
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TEST 2 



CIRCUIT DIAGRAMS AND DATA 
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B. 



TESTS 3-6 



The fixed exponentiation of a variable was evaluated 



in Tests 3-6. Tests 3, 4, and 5 were failures as the 
essential concept of current division to achieve the 
exponentiation was not understood. Test 6 demonstrates the 
use of a resistive network which accomplished the current 
manipulation. It is now apparent that as a result of the 
characteristics of the simple log configuration, current 
division can be easily accomplished by assigning the 
appropriate resistor value to the antilog input. This 
is diagrammed: 
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TEST 3 

CIRCUIT DIAGRAMS AND DATA 
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TEST 4 
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TEST 5 



- CIRCUIT DIAGRAMS AND DATA 
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TEST 6 

CIRCUIT DIAGRAMS AND DATA 
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C. 



TEST 7 



It is desired to validate the circuit segment for the 
operation 

x - ( PB C ~- c ' 



ec 



The operational sequence is as follows: 





The circuit described in the test data was set up to 
evaluate the feasibility of the operation. 

The desired operating characteristic of the cascaded 
logger - antilogger is picutred as follows : 
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The first log geneator establishes a voltage - 
current relation based on the diode and 20K resistor. 
Consider the second log generator to be superimposed in 
the first with a 3K resistor and a second diode. 

The cascaded antilog - configured amplifiers 
displayed gross thermal sensitivity. The simple OP amp - 
transdiode configuration is unsuitable for this cascaded 
operation without additional temerature compensation. 
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TEST 7 



CIRCUIT DIAGRAM 
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D. TEST 8 

i . Purpose 

Test 8 was conducted to verify the variable exponent 
part of the analog circuitry, using the four quadrant 
multiplier. 

The equation of interest is written below: 



CC 
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ii. Methodology 

The functional method used to exponentiate to a 
variable power was to regulate the multiplier output 
current in the ratio called for by the exponentiation. For 
example : 
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For multiplier scaling the integral buffer OP amp 
in the 2208 was used for the current division with a 
resistive bridge as follows, 




such that i-^ = i 2 and the output voltage of the multiplier 
is thus controlling the current input to the antilog 
amplifier transdiode. 

The multiplier output is nominally: 
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It was determined that a resistive bridge must be 
used for the multiplier input from logging OA, sketched 




as follows: 
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so by modulation of CC input, the desired current ratios 
across multiplier can be achieved. 

iii . Data 

A 5 minute observation of the steady state output 
voltage was made to investigate the circuit thermal 
stability. 
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10 Second Samples -Volts x 10 
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E. 



TEST 9 



i . Purpose 

To investigate the Texas Instrument TL 441 log 
amplifier. 

ii . Log Configuration 
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The output characteristic ^ L-oC» Avr v tc>w A>z 

verified. From the plot, "TL 441 transfer characteristics, 
the log amplifier has 1 1/2 decades of dynamic input range, 
from 20 - 600 mV. 

iii. Antilog Configuration 
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Antilog Configuration 1 
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The antilog data is graphed in "TL 441 antilog 
characteristics" for configuration 1. Configuration 2 will 
expand the antilog dynamic range somewhat above the 1 1/2 
decades available for configuration 1. 
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F. TEST 10 

i. Purpose 

To incorporate the 441 log amplifier and the 
quadrant multiplier in the variable exponent circuit. The 
general wiring is sketched below: 





INN 



ii. Steady State Performance 

The circuit demonstrated the ability to track a 
desired variable exponent. It was of primary interest 
to determine if the thermal properties were improved over 
the simple log-fed multiplier circuit. 
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The following data was generated: 
Variable Exponent Test 
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This represents a 50% reduction in the drift over 
the simple log configuration. Since the 441 log amplifier 
is fairly cheap, it is substituted for the simple log 
configuration. 



- 93 - 



APPENDIX IV 
DYSYS SIMULATION 
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APPENDIX V 

DESIGN DRAWINGS AND SKETCHES 
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APPENDIX VI 
SPECIFICATIONS 



A. 



ELECTRONIC COMPONENTS 



1 . Operational Amplifiers 
a. Fairchild X.A 741 



Maximum ratings 
supply voltage 

Power dissipation 

Supply current 

Typical input offset 
voltage 

b. Fairchild /°A 324 Quad OP amp 

Maximum ratings 
supply voltage 
single polarity 

Input voltage range 

Power dissipation 

Supply current 

Typical input offset 
voltage 

2 . EXAR 2208 Four Quadrant Multiplier 

Maximum ratings 
supply voltage 

Power dissipation 

Typical supply current 

Typical output nonlinearity 



+ 18v 
500 MW 
2 . 8 MA 

1 mv 



32v 

-0 . 3v to 32v 
900 MW 
2 MA 

2 mv 



+ 18v 
750 MW 
5 MA 
. 5% 
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3 . Texas Instruments TL 441 Logarithmic 
Amplifier 



Maximum ratings 
supply voltage 


+ 8v 


Input voltage 


6v 


Power dissipation 


500 MW 


Typical differential 
output offset voltage 


4 0 Mv 


Typical supply current vcc + 


18.5 MA 


vcc 


8. 5 MA 



4. Trimming Potentiometers 

All applications: Spectrol type 64 

ceramic potentiometers, l-100Kfl 

5 . Resistors 

All applications: Corning 2% 

1/4 watt resistors 

6 . Power Supply 

Analog Devices model 975 dual polarity 
power supply 

Output voltage + 15v 

Maximum rating 500 MA 

Typical output error 1% 

B. PERIPHERALS 

1. Breadboard Panels 



Continential Specialties model QT 59 
6.5" panels 
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2 . Volt Meters 

Beede Instrument 3" panel meters, 
model 3-03-8 

o 

250 movement 

200ft internal resistance 



3 . Joystick 

2 - axis movement; two, lOOKft rotary 
potentiometers on each axis. 



4. Chassis 



.062" clear anodized aluminum assembled 
with pop-rivet construction. 
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